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Abstract 


A theory of electromagnetic induction within a hemi-spherical conducting sheet 
is described. The theory is applied to the induction by Sg in a large ocean of 
uniform depth bounded by two meridians. The patterns of the induced currents 
are obtained and illustrated both for the 24- and 12-hourly components. Unlike 
the previous theories, the effect of self-induction, which turns out to be rather 
important, is fully taken into account. The magnetic fields produced by the 
currents induced in the ocean amount to several gammas. The electric field in 
the sea associated with the induction would be of the order of mv/km. 


1. Introduction 


The influence of seas and oceans on transient geomagnetic variations has been a 
topic of geomagnetic study for many years. Since it is highly probable that, in some 
cases, fairly intense electric currents are induced in sea-water which has an electrical 
conductivity of the order of 10°"'e.m.u., the influence may affect the interpretation of 
the relationship between both the parts of magnetic field originating from outside and 
inside the earth. 

S. Chapman and T.T. Whitehead (1923) made a mathematical study on the possible 
influence of a hypothetical ocean covering the whole surface of the earth with a 
uniform depth. They found that the internal part of a geomagnetic variation is con- 
siderably affected even by a shallow ocean. It seems unlikely, however, that the in- 
fluence of the real ocean is as large as that studied by them because the irregular 
distribution of land will greatly reduce the influence concerned. B.N. Lahiri and A.T. 
Price (1939) also made a similar sort of investigation with special application to the 
electromagnetic induction within a non-uniform earth. 

Since considerable influence of an ocean spread all over the earth has been found 
likely, it becomes important to estimate to what extent geomagnetic variations will be 
influenced by the actual distribution of land and sea though an accurate estimate would 
be of great difficulty to make. J.M. de Wet (1949) calculated the electric currents 
induced in oceans by Sg. On the basis of the distribution of ocean depth, he constructed 
a world-map of the distribution of the electrical conductivity over the earth’s surface, 
the boundaries between the lands and seas being approximated with suitable meridians 
and parallel circles of every 15°. Ignoring the effect of self-induction, he calculated 


( 65 ) 
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the induced electric currents by making use of the relaxation method getting current 
votices in large oceans such as the Pacific, Atlantic and Indian. 

T. Rikitake (1950) has also investigated the possible influence of a sea which is 
bounded by two meridians 72/2 apart. It was concluded that the electric currents 
induced by Sq would produce a magnetic field of several gammas. Since, the self-in- 
duction is completely neglected in these two studies, however, the results are not very 
accurate. A.A. Ashour (1950) studied the induction of electric currents in a uniform 
circular disk. He applied the result to the problem of the electric currents induced in 
an ocean which is regarded as a large disk. 

In relation to the anomaly of geomagnetic variations of short period in Japan, ake 
Rikitake and I. Yokoyama (1955) made a mathematical investigation of the induction 
of electric currents in a hemi-spherical conducting sheet placed in a time-dependent 
magnetic field as well as that in a plane sheet. These studies are later extended to 
the case of a spherical conducting sheet having a non-conducting hole (Rikitake eve 
1959). Although these studies have nothing to do with the electromagnetic induction 
by Sq, the importance of self-induction has been well demonstrated through the work. 

It has not been known to what extent the induction by Sg will be affected by self- 
induction. One of the purposes of this paper is to clarify such an effect. Assuming 
that a large ocean such as the Pacific is approximated by a hemi-spherical uniform 
sheet bounded by two meridians, a theory of electromagnetic induction by the 24- and 
12-hourly components of Sg is studied in this paper not ignoring the self-induction. 
The theory of electromagnetic induction within thin sheets or shells has been developed 
most rigorously by A.T. Price (1949). The theory that will appear in this paper more 
or less follows Price’s way of investigation. 


2. Inducing field of Sq 


It may be assumed, with a high degree of approximation, that the mean variation 
of Sq depends solely on local time, so that the magnetic potential of the inducing field 
arising from outside the earth can be written as 


W.= ay Di(r/a)"(en, ecosmat + en,s;sinmat)P"'(cos6,), allie) 


where a, t and 6) denote respectively the earth’s radius, local time and colatitude. a 


is the factor for the angular measure of the time and is given by aaa em 
60% 60x24 eos 


and é;,,; can be computed from the results of spherical harmonic analysis of Sg. P” is 
Schmidt’s spherical surface harmonic. 


If we denote the ¢.=0 meridian time by ¢, we have 
at=at+¢o. ( 2.) 
Putting (2) into (1), (1) can be written as 


W, =aX) X(7/a)" {en,ecosm(at + bo) + e@n,sSinm(at + ¢o)} P”(coso) (CB) 


which, by adopting complex expression, becomes 


. 
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W. = aX) di(7/a)"(Criy af ie”".Jet mCat+ #) P™(CcOsG) ; ( 4 ) 


We may hereafter take (4) as the inducing field of Sg under the condition that suitable 
account is taken of the real and imaginary parts of the final results. 


Zo In the above, the potential has been 
defined in reference to a coordinate of 
which the @)=0 or zp) axis is taken as the 
line connecting the centre of the earth to 

Y(Y) the north pole as is shown in Fig. 1 in 
which the hemi-spherical sheet which is 

Xo(z) supposed to be the ocean is also shown. 

In order to obtain solutions of good con- 

vergency, however, it seems better to 

transform the coordinate so as to let the 

pole agree with the centre of the ocean. 

As canbe seen in Fig. 1, the x0, Yo and 2 

X axes are transformed to the x, y and z 

Fig... 1 ones respectively. This sort of transfor- 

mation is useful for limiting the orders of 

spherical surface harmonics which will come out in the course of the calculation. 

Actually, only P,! and P,? will be needed for studying the induction by the 24-hourly 

component, and P,', P,? and P,? for the 12-hourly one as can be seen in the following. 

Y. Sat6 (1950) studied transformation of spherical coordinates in a general way. 
Applying his results to the present case, we obtain 


P,1(cos0y)cos¢ y= — P;'(cosd)cos¢, 
P,1(cos69) sing y= P;7(cosA) sin2¢, 


P,2(cosb)c0s26)= ~/10P,"(cosd)cosp-+ 45 P,%(cosd)cos3¢, (5) 


P22(cos09) sin2¢)= — P3*(cos@) sin2¢, 
where the quantities with subscript 0 refer to the original coordinates. 
It has been well known that harmonic terms such as P,' and P;? are most im- 


portant for the mean Sq, so that we may assume that the inducing field can be divided 
into two parts only neglecting harmonic terms of higher order. We write 


W.=W,+ W;,, (6) 
where 
W,=a(r/a)*(éh,¢ —163,,)@* "+ © Po (COSA), | (7) 
W,=a(r/a)*(G,0 — i€3,)e?* °° P3*(COSd,) « 


Applying the transformation formulas which are given in (5) to(7), we get 


W,=a(r/a)’e*™(@,¢— ieé},;)(— P2'(cos0)cos¢ + iP.2(cos)sin2¢), 
We=ale/aye***(G,6— id.)( 10 p,s(cosd)cosp —iP;?(cosé)sin2¢+ 6 P,*(cos6)cos36) , 
(8:) 
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which can be regarded as the inducing field of Sq. @i,c, Cassy é,, and ¢@,, are to be 
obtained referring to the ¢)>=0 meridian being the time origin. If we suppose that the 
Pacific Ocean is approximated by the hemi-spherical sheet and identify the ¢o= 

meridian with the 120° E meridian which approximately agrees with the west boundary 


of the Pacific Ocean, these coefficients are easily calculated from the analysis of Sg 
If we take Chapman’s values (Chapman S., 1919), they become as shown in Table 1. 


Table 1. The coefficients of the harmonics of magnetic potential of Sg in units 
of gamma referring to the 120°E meridian being the time origin. 


| C555 


| —5.8 


| aye 


| 3.2 


1 
C 208 


Cs,¢ 


5.6 | 5.0 


3. Current functions of the induced electric currents 


Corresponding to the inducing fields W, and W,, the current functions of the electric 
currents induced in the hemi-spherical sheet, which occupies the part of the surface 


A A a 


from 0=0 to 0=7/2, can be written as 


V =e! +{ Ya,P2(cosd)cos¢ + 310,P2(cosf) sin2p | 
(9) 


v, ag! “| SanP2(cosf)coss + L1BnP2(cosd)sin2g +3 rn *(cost)cos3¢ | 1 


It should be remarked that these expressions hold good not only on the hemi spherical 


‘sheet but also on another half of the surface (z/2<@<z). 
According to the theory of electromagnetic induction within thin sheets (Price A.T., 


1949, Rikitake T. & Yokoyama I., 1955), the scape that should be satisfied by the 


EIELER uae on the conductor is at 


ec 2 _ bor Baie) #br aan Ste ele , : ge a fs 
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where 


ees (13) 
2raka : 


4. Induction by the 24-hourly component of Sq 


Multiplying by Py! sin@, we integrate (11) with respect to @ from 0=0 to 0=n, 


whereas the first expression is taken for 0<@<z/2 and the second one for 2/2<0<z. 
We then obtain 


—MN+1)ayRyy=@Kiolei,, —iti,s)Roy+ ig AE aan (14) 
n 2n 
where 


Rae Wesnea'y “(ede . (15) 


Rnry for nZ6 and NZ6 are given in Table 2. 


Lable:2)) Ran (=Ryn) 


a (14) are the simultaneous equations for a,’s solving which we can determine them. 
In a similar way, | we multiply (12) by Py sin@ and intectals its fre simultaneous = 
~ equations for bx'8 are then written as BATE ee oH S24): % oe — 2 
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If we denote the real and imaginary parts of @, and 5, respectively by a, and a* 
and b, and 6,*, (14) and (16) lead to 


SENNA Din Rey mt 1) a eee 


—EN(N+lay*R SAS Spe a) a 
ENN-+ 1)ay* Ryn — cag a on 2N 
and 
~§ NN+1)bySww $I MOED 6,889 = 8h* Soy, 
n-- 1 2 | 
n(n+1) antes. saa | (19) 
—EN.N+1)by*Syv— Pa eee caer = Sant 


In the case of Sg, € amounts to 8.59, 0.859 and 0.172 respectively for the ocean of 
which depth is assumed to be 100, 1000 and 5000 m, while the conductivity of sea-water 
is taken as 4X107!'e.m.u. Adopting 1000m as the depth, the simultaneous equations 
(18) and (19) are obtained as given in Tables 4.and 5 where harmonics for ~>6 are 
ignored. As to e},, and e},,, the values given in Table 1 are used. 

Since the diagonal terms of these simultaneous equations are fairly large, they 
can be easily solved. The solutions in units of gamma are given in Table 6. The 
figures in the table should be multiplied by @/2z. 


Table 6. The coefficients for the induction by the 24-hourly component of Sq. 


ee dt 
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a,* | a | a,* | a,* | a,* | ag* 
2.118 | 1.087 | 0.336 | 0 | —0.042 | 0 
a by ays ONS or vines 
1.087 | 0.553 0 —0.109 | 0 
bo | byt by bX | bX 
i 0.970 | _0.475 0 0.098 | 0 


; Comparing Table 6 to Table 7, considerable differences are found between the cor- 
responding coefficients. We therefore see that the effect of self-induction is so large 
that no accurate results can be obtained when it is ignored. If the conductivity of 
sea-water or the depth of the ocean are taken to be larger than those of the present 
case, the influence of self-induction becomes larger. 

The coefficients given in Table 6 enable us to obtain the current function and 
consequently the magnetic potential or field produced by the induced currents. The 
condition that the component of the electric currents normal to the ocean boundary 
should vanish at 6=z/2 leads to an alternative expression of the current function being 
written as 


Ua «| DcamPsn(cosd)cosg + didom- iP ameCOSE) sin26| <e (20) 


which is true only for 0<@<z/2. On the other half of the sphere, ¥; is zero everywhere. 
We then have tae ase = 


— 


F- SasP i= TeenPie 00K 0/2, | at 
DesPi=0 7/20 <i ees) ae ei = 
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are obtained from those in Table 6, the electric current systems induced in the hemi- 
spherical ocean are easily obtained on the basis of (20) for epochs at=0°, 45°, 90° and 
135° as are respectively shown in Figs. 2, 3, 4 and 5. In these drawings, the current 
systems are viewed from a distant point right above the centre of the ocean. Electric 
currents amounting to 2500 amperes are flowing between adjacent lines. These epochs 
correspond to 0, 6, 12 and 18h of the ¢)=0 meridian time. 

The magnetic field produced by the induced currents can be also estimated by use 
of the coefficients given in Table 6. For illustrations, the northward components at 


——----— 


_—— 


Seles, 


SP miles 2 “SP 


Fig. 2 Induced currents for af=0° in the hemi- Bigs 3) ai=45° 
spherical ocean by the 24-hourly component of Sq. 
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Fig. 6 The northward component of the induced Fig. 7 The vertical component of the inducing 
magnetic field at at=0° along the meridian pas- magnetic field at at=0° along the meridian pas- 
sing respectively the central points of the ocean sing respectively the central points of the ocean 

(upper curve) and the continent (lower curve). (upper curve) and the continent (lower curve). 


at=0° are calculated along the meridians connecting the poles to the central points of 
the ocean and continent. They are shown in Figs. 6 and 7 from which we see that 
the magnetic field produced by the currents induced in the ocean amounts to about 6 
gammas at the maximum which occurs on the equator for the northward component 
and at both the centres of current vortices for the vertical one. Meanwhile, the magnetic 
fields on the continental or non-conducting side of the sphere are so small that they 
are probably within the range of errors of estimation. The eastward component of 
magnetic field becomes large near the boundary of the ocean, its maximum magnitude 
being also several gammas. We thus see that the effect of the electric currents induced 
by the 24-hourly component of Sq in a hemi-spherical ocean 1000 m in depth is some- 


thing of several gammas. 


5. Induction by the 12-hourly component of Sq 


The induction by the 12-hourly component of Sq can be studied in the same way. 
The simultaneous equations for the real and imaginary parts of a,’S, Bs’. and 7,’s, by 
which we may calculate ¥, by the second equation of (9), can be obtained as 


—ENN+ l)oyRyy +231 aU Don Noa Sang : 1065, Roy , 
ae (25) 
—ENN+1)ay*Ryy— Re ae pen Ran ar, 31064, Ray 
_EN(N+1)By Suey +2520) 8, Sany= —2 Bel, 0Ssw 5 
nm Qn+l1 2 (26) 


r= ( # =e lg EA) Brn = ieul 3€5,0Say » 
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HEN lity Tan eee oe MTD Tg Tye 
gant _ (27) 
—~EMN+Dry*Txw—-25 mi WW ieey perny ea? Si pe 
2n+1 1 
where 
ee | Pl) PHa)de (28) 


T,y for nZ6 and N46 are given in Table 8. 


Table 8 Tay(=Twn) 


as 3 4 5 6 
3 0.2857 0.2057 0 —0.0716 
4 : 0.2222 0.1122 0 
5 0.1818 0.1184 
6 0.1538 


Solving (25), (26) and (27) in a way similar to that in the last section, the coefficients 
are determined as are given in Table 9. The values for e,. and é,, given in Table 1 
are used in the numerical calculation. The units of the values that appear in Table 9 
are the same as those in Table 6. 


Table 9 The coefficients for the induction by the 12-hourly component of Sq. 


a a 
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Fig. 8 Induced currents for af=0° in the hemi- Fig. 9 at=22.9° 
spherical ocean by the 12-hourly component of Sq. ; 


NP. NP 


ean SP 
Fig. 10 at=45° Bigs 10 ai—67.97 


larger than that for the 24-hourly component. We would have considerably different 


patterns of induced currents if the effect were ignored. The coefficients for the case of 


no self-induction are given in Table 10. 


Table 10 The coefficient for the induction by the 12-hourly component of Sq. 
The self-induction is ignored. 


ay | Ao | a3 | a4 | as | a6 | a;* ast | ast at | ast | ag* 


EO | —1.180| -1.334] —0.654 | 0 | 0.145 | 0 | 9-651 | 0.736 | 0.361 | 0 | -9.080) 
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| Bo a) Bs | ; Bs | Bs Be a B3* jee B* joes 

| ‘ ee a eee eee ea | Raa ll 5 

1.302 0.931 | 0.340 | 0 |-0.050] 2.360] 1.668] 0.616 | 0 |—0.091 
| 73 | 7 | if Te rs* pe veie ree 
| 1.034 | 0.579 | 0 | —0.187| —0.570| -0.320| © 0 0.076 


6. Discussion and conclusion 


A theory of electromagnetic induction in an idealized ocean is developed in this 
paper. Applying the theory to the induction by Sg ina hemi-spherical ocean of uniform 
depth bounded by two meridians, the patterns of induced currents are obtained both for 
the 24- and 12-hourly components of Sq. 

It is found that the patterns change according to the universal time and the magnetic 
fields produced by the currents amount to several gammas. From the illustrations 
showing the current patterns, we see that the maximum value of the induced currents 
is roughly 10000 amperes per 10° in latitude, so that electric current of the order of 
10-4 ampere/cm can be induced by Sq in the ocean 1000 m in depth. The current density 
in that case becomes 107° ampere/cm?. Hence the electric field in the sea is estimated 
at 2.51078 volt/em because the conductivity has been taken as 4x107!!'e.m.u. If we 
suppose that two electrodes are put in the sea, the distance between them being 1 km, 
the potential difference between them amount to 2.5mV. For an order-of-magnitude 
estimate, therefore, the voltage induced in an ocean by Sg would be a few millivolts 
per km or thereabout. 

Although no account has been so far taken of the effect of self-induction in studying 
the electromagnetic induction in an ocean by Sq, it turns out that the effect seriously 
affects the induced currents. The current functions obtained without taking self-induction 
into account differ considerably from those with proper account of it. It is emphasized 
that the effect of self-induction should be fully taken into account in studies of this 
sort. 

Since the magnetic field produced by the electric currents induced in a large ocean 
amounts to about 1/3 of the inducing field, the Sg seems to be affected by such an 
ocean to some extent. It is not known, however, how and to what extent the actual 
oceans will affect Sq because of the complicated distribution of lands and seas. It will 
be highly desirable to detect influences of large oceans on Sq by conducting a close 
analysis of observed data from well-distributed observatories. This sort of investigation 


would be also necessary for accurate determination of the overhead current systems of 
Sq. 
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Physical State of Outer Atmosphere and the Origin 
of Radiation Belts* 


By Tatsuzo OBAYASHI 


Radio Research Laboratories, Kokubunji, Tokyo 
(Read November 27, 1959; Received January 22, 1960) 


Abstract 


A possible mechanism trapping high energy particles in the radiation belts 
surrounding the earth is proposed, taking into account the existing hydromagnetic 
waves in the outer atmosphere. It is shown that there are two regions where 
the amplitude of hydromagnetic waves and the compressibility of gas bearing 
magnetic field are large. Since the acceleration mechanism is operative in such 
regions, particles inside the ‘regions may be raised in their energies, and con- 
sequently they produce local inhomogeneities of hot plasmas. These hot plasmas 
interact with the geomagnetic field and may form a certain kind of magnetic | 
bottles, in which the high energy particles are likely trapped. The inner activated | 
region is expected at the height 1,000~3,000 km and the outer region is of the order : 
of 20,000km. Although the origin of high energy particles is possibly of injected 
solar particles, their concentration at particular regions may be controlled con- 7 
siderably by the hydromagnetic nature in the outer atmosphere. 


1. Introduction 


During the past few years, considerable advance have been achieved in the 

knowledge of the upper atmosphere beyond the ionosphere. New evidence indicating 

the existence of regions trapping high energy particles has been discovered by Van 
2 Allen et al [1958], which are now called the radiation belts. These are two distinct 
= ; Be ae Ot see tc: verticeorithesinnes, caneslencth ast a fore tots ; a 


er 
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appear due to the interaction of fluid motions with the magnetic field. Waves, 
known as Alfvén waves, are propagated along the magnetic lines of force with a 
velocity. 
é H 
AS. 7 
/ 410 
whepe H is the magnetic field intensity and o is the ionic density in the medium. 


Besides this Alfvén wave, hydromagnetic waves have two other modes; the modified 


C1) 


Alfvén waves and the retarded sound wave. In particular, for the propagation of 
waves parpendicular to the magnetic lines of force, the velocity of the Alfvén wave 


and the retarded sound wave is zero, but the modified Alfvén wave has a velocity 
v=,{ #4 Ve (2) 
Arto 


where V,=./2kT/m, sound wave velocity. This wave, in which the particle velocity 
is parallel to the direction of propagation, both being perpendicular to jg, 1S 2 


longitudinal hydromagnetic wave. 

There is a possibility that such waves may be generated in the earth’s outer 
atmosphere due to excitation by solar corpuscular streams, and produce hydromagnetic 
oscillations. It has been shown elsewhere [Obayashi, 1958] that there exist two kinds 
of mode of oscillations in the outer atmosphere. The one is caused by a transversal 
hydromagnetic wave (Alfvén wave) propagating along the geomagnetic line of force, 
whose oscillation is perpendicular both to the direction of wave propagation and the 
magnetic field, and each surface of revolution of a geomagnetic line: of force oscillates 
independently. The other is the longitudinal oscillation, the disturbances propagate 
across the geomagnetic lines of force as the longitudinal wave mode, and the oscill- 
ation is set up between the top of the ionosphere and the outer boundary of the earth’s 
atmosphere where solar plasmas are invading. 

There is considerable evidence to believe that geomagnetic pulsations observed at 
high latitudes are mainly the transversal type, while at low latitude the oscillation 
of longitudinal waves prevails. The longitudinal hydromagnetic waves may have 
an important effect on the formation of radiation belts, because of their heating or 
accelerating actions of atmospheric gas. An idea, that hydromagnetic waves are res- 
ponsible for producing high intensity particle radiation, has been first suggested by 
Dessler [1958], and he has shown that existing large hydromagnetic waves with sharp 
crests above 1000km are favourable for the Fermi type particle acceleration process. 
However, as will be described, it is shown that this mechanism producing high energy 
particles is not sufficient enough to account for the observed energy spectrum of 
radiation belts. Therefore, some re-consideration of this theory will be made in the 


following section. 
3. Origin of Radiation Belts 


In an agitated plasma bearing magnetic field, high-speed charged particles gain 
their energies by colliding with moving clumps of magnetic field or hydromagnetic 
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waves, and thus undergo Fermi acceleration [Fermi, 1954]. According to Parker [1958], 
the effectiveness of the Fermi acceleration depends on two factors, one is the amplitude 
of hydromagnetic waves and the other is the compressibility of the plasma. As Has 
been pointed out [Obayashi, 1958], there is an evidence that longitudinal hydromagnetic 
oscillations exist in the outer atmosphere which are excited by solar corpuscular streams. 
The mode of such oscillations depends on the ionic density and the geomagnetic field 
distribution in the outer atmosphere. Using the current model of the outer atmosphere 
shown in Fig. 1 (due to rockets and whistler observations), the velocity of a long- 
itudinal hydromagnetic wave is 
computed using eq. (2). It is clear 
that there are two regions where 


1000 5000 10000 20000 30000 Height (Km) 


5 the wave velocity decrease with 

S Electron Density Distributio i ; 

2 height. As has been mentioned by 

= y Whistlers Dessler, in these regions downward 

s | F 

. najoten TR) il eeio travelling hydromagnetic waves 

> 7 will be reflected back. Hence, it 

; Ho=100Km nal Ss ‘ 
r= seioe y=24) NoseWhistlers is likely that the hydromagnetic 
oscillation may attain a large ampli- 

= tude at the fringe of the ionosphere 

3 pices i 

8 “; near the height of 1,000km and 

$ Modified Alfven Wave 


| above 10,000km while it has a 
2 V4 C 

S=(Se/nkT)->104_ (Mag.Mach Number) kind of nodal point in the ionos- 
Effective Particles Acceleration phere and at 5,000 km above the 


He Apc a earth’s surface. 


On the other hand, the Fermi 
process can be operative only as 


long as the hydromagnetic waves 
can maintain sharp fronts or crests. 
This situation 1s possible only when 


the compressibility of the gas bear- 


a Ascending Flight . pin , 
\y Descending Flight ing magnetic inhomogeneities[ 4H? 


/8xnkT] is large. In Fig. 1, the 
region of the magnetic Mach 


Pioneer IL 
Van Allen Bands 


counts /sec. 


7. number larger than 10? is indicated 


Distance in Earth Radii by the line. Thus, it can be under- 
Fig. 1 Hydromagnetic oscillations in the outer 


ae stood that th i 
atmosphere and the position of the radiation belts : pe AIS Lat SI 


where the situation is favorable 

for atmospheric gas particles to be accelerated by hydromagnetic waves. As shown 

by hatching in Fig. 1, the region I exists at the height between 1,000~3,000 km, and 

the region II is of the height between 10 ,000~20,000 km. It is worthwhile to note that 
they are situated at the place where the radiation belts have been found. 

The energy spectrum of accelerated particles by the Fermi mechanism (non rela- 
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tivistic) is given [Parker 1957], 


aN il w-Ww 
j(W)=- 3 |- = 2] 
Mf) 2N mV? aR 2NomV? - oe 


where W, is the energy of thermal ions in the outer atmosphere. N, is number of 
collisions of a particle with hydromagnetic wave of the wave length 4 before escaping 
from the agitated region of dimension L. If particles undergo random walk in the 
region, Nj is given by (L/4)?. Assuming that V=1,000 km/s, 4=5,000km andZ=10,000 
km it can be shown that the energy of protons attained by this process is 


W*~2NmV?~100KeV ‘si 


Therefore, the Fermi acceleration is not sufficient to produce the observed energy, although 
if an equipartition of energy between the electrons and the accelerated protons is 
achieved, the existing electrons in the outer radiation zone may be explainable by this 
mechanism [Singer, 1959] 

However, since the regions activated by hydromagnetic waves are situated at the 
place where the radiation belts have been found (see bottom of Fig. 1), the following 
discussion may be appropriate: Owing to the continual agitation by hydromagnetic 
waves, part of the particles inside the activated regions may be raised their energies. 
Those high energy particles will suffer collisions and change their kinetic energy to 
heat. The heating may also be due to the violent beating of neutral gas by hydro- 
magnetic waves. Consequently, the continual heating of the outer atomsphere by such 
processes may generate local inhomogeneities of hot plasmas. It is suggested that 
those hot plasmas expand outwards interacting with the existing geomagnetic field. 
The geomagnetic configuration may be altered to some extent like a magnetic bottle, 
shown in Fig. 2. Two regions of such magnetic bottles are expected, where high 
energy particle are trapped more firmly. This 
may be consistent with the result of the 
measured outer geomagnetic field by the 
“Cosmic Rocket” [Dolginov and Pushkov, 
1959], which shows a slight drop of scalar 
field intensity below the theoretically expected 
one in the radial distance range 15,000~22,000 
km. 

The origin of high energy particles is 
still uncertain. However, observation by 
Pioneer IV, which was launched after a geo- 
magnetic storm, indicates an increase of 


particle flux even at the inner radiation belt 
[Van Allen and Frank, 1959]. The height 
fii meee eee oe of the inner belt decreases considerably during 
magnetic storms [Yoshida, 1959]. These 


ne suggest that the injection of solar particles is likely the origin of radiation 
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belts, and also that the activation by hydromagnetic waves may take an important 


role to maintain and control the radiation belts. 
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Magnetic Susceptibility of Some Rock Forming Silicate 
Minerals Such as Amphiboles, Biotites, 
Cordierites and Garnets* 
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. Abstract 


Some silicate minerals such as amphiboles, biotites, cordierites and garnets 
were examined magnetically. Magnetic susceptibility at the room temperature 
of garnets is in good agreement with the theoretical value estimated by using the 
effective Bohr magneton number derived empirically from salts bearing above 
magnetic ions, while in the case of amphiboles, biotites and cordierites, which 
contain water in their crystal structure, magnetic susceptibility at the room tem- 
perature shows fairly larger value than expected. 


_— 


1. Introduction 


Magnetic properties of ferromagnetic_ minerals contained in rocks have been ex- 
amined fairly comprehensively, but only few of the silicate minerals such as olivines 

_ (Nagata, Yukutake and d Uyeda : 1957), orthopyroxenes (Akimoto, Hérai and Boku: 1958), 
and monoclinic pyroxenes (Chevallier et Mathieu: 1958 and 1959) have been examined 


WD 
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Table I Chemical analyses of amphiboles, biotites, cordierites and garnets 
a) Amphibole 


| No. 1 No. 2 No. 3 No. 4 No.5 No. 6 
SO, “| Swi 43.20 44.07 54.09 56.61 45 
ALO, | 0.61 12.44 12.37- ~ 0.22 0.24 4 
TO, | 0.45 1.65 1.70 3.54 1.04 0 
Fe,0, | 14.37 3.21 0.18 0.66 1.10 10 
FeO ) 14.82 10.10 10.23 4.69 3.90 26 
MnO | 1.46 0.21 0.18 0.26 0.29 0 
MgO 5.07 13.27 14.20 21.29 22.27 0 
Go “Pris 11.36 12.42 12.14 12.33 2 

| NaO | 4.94 2.72 1.00 0.21 0.06 7 
K,O 2.10 0.40 0.30 0.16 0.18 3 
H.O* | 2.02 1.38 2.85 0.12 0.10 1 
H,O- 0.10 0.07 0.03 2.78 2.45 0 


99.53 100.16 100.57 


Mineralogical description 

No. 1 Riebeckite from quartz-syenite, (AM 440702-10), Loc. 500M. W. of Bungan, 
Koganzan, Heiko-yu, Heiko-gun, Kogen-dd, Korea. (Analyst ; Inoue) Ref. 1 

No. 2 Hornblende from biotite-amphibolite, (AM470808-3), Kamata, Hurudono- 
mura, Higasi-sirakawa-gun, Hukusima-ken, Japan. (Analyst ; Haramura) Ref. 2 

No. 3 Hornblende from epidote-hornblende-plagioclase-schist, (AM470621-16), 
Kamiyama, Tabito-mura, Iwaki-gun, Hukusima-ken, en grey sical 
shiro), Ref. 2 — ' ; bs 

No. 4 Acie dine Chesten Vermont, USA. Ref. 3 


SS ene i oe - oe 


"The ae of No, 8 ae rears ones acing 1 Dr anos bites 


Magnetic Susceptibility of Some Rock Forming Silicate Minerals 87 


Mineralogical description 

No. 1 Biotite from cancrinite-syenite-pegmatite, (AM 440617-14 (11)), Loc. Fuku- 
shin-zan, Sei-men, Heiko-gun, Kogen-d6, Korea. (Analyst ; Inoue) Ref. 1 

No. 2 Biotite from nepheline-syenite, (AM 431205-13 (7)). Loc. Shacho-zan, 
Kennai-men, Heiko-gun, Kogen-do, Korea. (Analyst ; Inoue) Ref. 1 

No. 3 Biotite, (AM 470623-7), Between Isizumi and Kaiya, Tabito-mura, Iwaki- 
gun, Hukusima-ken, Japan. (Analyst; Haramura) Ref. 2 

No. 4 Biotite, Specimen of Geological Institute, Tokyo University. 

No. 5 Biotite from nepheline-syenite, (AM 440619-13 (12)), Loc. Fukushin-zan, | 
Sei-men, Heikd-gun, Kogen-dd, Korea. (Analyst ; Inoue) Ref. 1 

No. 6 Biotite from garnet-biotite-oligoclase-quartz-gneiss, (AM 540812-6A), Nisi- 


takat6, Takaté-mati, Kami-ina-gun, Nagano-ken, Japan. (Analyst ; Haramura) 
Ref. 2 


Content of Fe,0;+FeO of No. 4 was determined from the ao sian index | 
by Dr. A. Miyashiro. 


—e— 


c) Cordierite 


Mineralogical description 
: - No. 1 Iron-rich cordierite from pegmatite, pura Oomorita- aT 4 
idle epee al Mladbe res ean ‘Hukusima-ken, Japan. (analyzed Hebi | 
eatin linitte icrevanatyekof Xt . ei 


it CEOS ah lira wae ca spent ae ea 
of ' ; 
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d) Garnet 

| No. 1 No. 2 No. 3 No. 4 No. 5 No. 6 No. 7 
SiO, | 39.68 36.82 38.53 38.74 

| Al,O, 22.05 22.10 veh © 22.09 

| TiO, 0.98 0.08 0.09 0.07 

| Fe,0, 4.10 1.54 1.66 2.97 : 
FeO 11.33 14.96 25.35 23.23 
MnO |g tamer 15.08 0.98 0.47 | 
MgO 13.63 2.73 8:12 14.20 
C20 6.75 5.93 3.21 1.10 
Na;O 0.19 0.12 0.18 0.23 . 

| K,O 0.04 0.09 0.04 0.11 | 

| H,O* 0.82 0.38 0.05 | 

; 0.28 

H,O- ) 0.00 0.06 0.02 . 

| P.O; | 0.01 

| Cr,0, | 0.70 

| Total | 100.60 100.28 100.72 100.35 


| 

Mineralogical description 

No. 1 Garnet from eclogite, Gongen-yama, esha Uma-gun, Ehime-ken, 
Japan. Ref. 8 

No. 2 Almandine, Ennis, Montana, USA. 

No. 3 Garnet porphyroblast in metagabbro, Gore Mt. (Barton mine), Adirondack, 
N_ Y., USA. 

No. 4 Garnet in garnet-quartz-nodule in biotite-schist, ae ee (Segregation 

nodule ?), Tonaru, Bessi, Ehime-ken, (Analyst ; Haramura) 
No. 5 Garnet from garnet-amphibolite (A-9), Prince Harold Coast, Antarctica, (An- 


alyst; Katsura) 

No. 6 Garnet from banded gneiss, (A-20), Prince Harold Coast, Antarctica, (Analyst ; 
Katsura) 

No. 7 Specimen of the Gebiogiat Jnstithibes Tokyo ieidehaits: 
There is no analysis of No. 2, 3 and 7 of garnets — a é 


Most samples" were originaly separated out by an isodynamic separator. However, — 
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where N, , and k are Avogadro number, Bohr magneton and Boltzmann’s constant 
respectively ; x, y, Z, and Pgr, Pay, Paz, are the number of the gram-ions of Fe**, Fe, 
Mn?*, contained in unit gram of specimens and its effective Bohr magneton number 
respectively. After Stoner (Stoner: 1934), the values of effective Bohr magneton number 
are taken as follows; Pyr=5.25 to 5.53 and pyy=Ppg.=5.88. x and 22 corresponds to 
Ppr=5.25 and 5.53 respectively. 

(b) Measurements at low temperatures 

Measurements at low temperatures were carried out by immersing specimens to 
be measured into oxygen liquid in the dewer. Some o-H curves are shown in Fig. 2. 

For the purpose of examining deviation of the magnetization curve from Curie’s 
law to Curie-Weiss’ law, the paramagnetic Curie temperature 9, given by x(T—@,))=C 
was obtained from the observed values of x at the room temperature and at low tem- 
perature, being also tabulated in Table II. 

(c) Variation of susceptibility with temperature 

In order to ascertain whether Curie’s law is valid or not, variation of magnetic 


Table II Magnetic susceptibility of amphiboles, biotites, cordierites and garnets, 
comparing with the theoretically calculated values and paramagnetic Curie 
temperature (unit; emu/gr). 


observed value 


Specimen Iculated 
P eolculaicd valve at room temp. at 90° K 
| 10°%, ‘Ca coe 10% 6,°K- 
Amphibole-1 516k 7.5 0.013 
2 2.45 2.2 0.010 
sl 1.91 33 0.001 
4 | 1.02 1.3 0.000 
5 | 0.96 13 0.000 
64 6.6 


8.0 0.000 
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susceptibility with temperature was examined ; namely y-7 measurement for the range 
from liquid oxygen point to about 500°C. Because all the specimens except garnets 
examined here become unstable above about 500°C, making decomposition and changing 
into magnetite etc., no measurement above 500°C was carried out. 

Examples of results of measurements are given in Fig. 3, where it is shown that 
x is inversely proportional to 7, but generally the y-!-T curve does not cross 0°K 
point. 


z 4, Discussions 


Amphibole - 6 


Summerizing above results, x 
at the room temperature~ x, and x2 
diagram and 0,~¥xz, and yx, diagram 
are given in Fig. 4 and 5. 

In Fig, “4, if7is “seen that othe 
© heating observed values of x at the room 


© cooling 


(in arbitrary scale ) 
T 


temperature are not in good agree- 
ment with the estimated values of 


<00" 600" rx 7, and y=, except for case of garnets. 


o 200° 
Fig. 3 variation of magnetic susceptibility The measured values of magnetic 
with temperature nee : Sate 

susceptibility of amphiboles, biotites 

and cordierites are larger than 
calculated values x, aud x, by 20% 


~30% generally. It seems that this 


Biotite -3 


discrepancy exceecs possible range 
of experimental errors, which can 
be estimated to be less than 10% 
at the largest. 


papicis Comparing chemical composi- 


© cooling 


(in arbitrary scale} 


tion of these minerals, magnetic 
susceptibility of which cannot be 


m4 me or coo vi exactly given by the Langevin’s law, 

with those of the minerals for which 

magnetic susceptibility is well represented by the Langevin’s law such as garnets, 

olivines and pyroxenes, it may be found that the former group minerals contain OH 

radicals or H,O molecules in their crystal structure, but in the latter there is no water 

component. In Fig. 5 also, it is shown that paramagnetic Curie temperature of the 
former group of minerals are fairly larger than that of the latter. 

This fact may suggest that some very weak interaction or something like weak 
magnetic coupling among magnetic ions through OH radicals or H,O molecules in these 
minerals containing water may affect their magnetic susceptibility. Further careful 
study: will be necessary in order to examine existence of such a weak interaction 


through OH radicals or H,O molecules, 
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Fig. 4 Relation between measured magnetic Fig. 5 Relation between paramagnetic Curie tem- 
susceptibility, y, at the room temperature and perature, 9», and theoretically calculated magnetic 
theoretically : calculated values, y, and x, , susceptibility, x, and x5 
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Abstract 


Physical mechanism causing the scattered polarization or the unstable 
magnetization of some dolerite sheets was investigated through various magnetic 
measurements and microscopic observation. It was found from X-ray and chemical 
analyses that titanomaghemite exclusively occurs in the rocks showing unstable 
magnetization. This implies that the scattered polarization or the unstable 
magnetization of the rocks may be attributable to the following mechanism: the 
thermoremanent magnetization is destroyed by natural oxidation from original 
titamomagnetite to titanomaghemite with the result that the remanent magneti- 
zation measured in the laboratory is only the isothermal remanent magnetization Z 
produced by the magnetic after effect under the geomagnetic field. The situations 
that the magnetic instability of rocks varies systematically according to the degree 
of the alteration or weathering were also elucidated. 


1. Introduction ; ~< 7 
_Tt has often been reported that there exist some rocks showing scattered magnetic 
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for the scattered polarization or the unstable remanent magnetization. Even the most 
fundamental question whether or not there is any difference between the ferromagnetic 
minerals in the rocks with stable remanent magnetization and those in the rocks with 
unstable one has been pending for years. In order to elucidate the physical origin of 
these phenomena, we have performed a study with a few magnetically unstable dolerite 
sheets. These dolerites were originally sampled for the palaeomagnetic studies of the 
Japanese Tertiary but found to be unsuitable for the purpose on account of their 
unstable magnetization. (Nagata et al., 1959) Common features to all these unstable 
rock specimens were carefully examined through various stability tests. Microscopic 
observation on the polished surface and X-ray or chemical analysis of the separated 
ferromagnetic minerals were also carried out. The systematic study along this line 
may in turn furnish us a more precise knowledge on the natural remanent magneti- 


zation of rocks and allow us to interpret the palaeomagnetic results more fully. 


2. Unstable Remanent Magnetization 


In the course of the development of our palaeomagnetic research on Japanese 
Tertiary igneous rocks, there have been found many rock masses having scattered 
polarization. Typical examples are the Miocene dolerite sheet swarm near Kingosi, 
Sidara basin, Aiti prefecture (sample nos. 2502 and 2503) and the Miocene dolerite sheet 
of Warabidaira in the valley of the river Huzi (sample no. 2601). In this section the 
characteristics of the unstable remanent magnetization of these specimens will be 
revealed through various stability tests. 

2. 1. Scattered polarization 

Six oriented samples were collected from a single dolerite sheet near Kingosi (2502) 
over a distance of about fifteen metres and from these fourteen specimen cores were 
cut out. The directions of the natural remanent magnetization of the specimen cores 
are shown in Fig.l-a with the aid of Schmidt projection. An appreciable scatter of 
the directions is found even in the specimen cores cut from the same rock sample. 

From another large dolerite sheet near Kingosi (about one kilometre apart from 
the locality of 2502) six oriented samples were collected over a distance of thirty metres 
(2503). Thirteen specimen cores were cut from these samples. The magnetic polari- 
zation for each core is plotted in Fig. 1-b. The widely scattered polarization was 
found in this site, too. 

For the dolerite sheet of Warabidaira (2601) which is about a hundred kilometres 
apart from Sidara basin, ten specimen cores were cut from six rock samples. Scattered 
direction of the remanent magnetization is shown in Fig. 1-c. 

In contrast to these scattered polarization, the directions as well as the intensity 
of remanent magnetization have been proved to show good consistency in many other 
rocks for our palaeomagnetic research of the Tertiary age. As an example, the 
directions of the natural remanent magnetization of the Miocene andesite lava flow at 
Tengudana, Sidara basin (sample no. 2401) are given in Fig. 1-d. This specimen is used 
in the present paper as the standard of the magnetically stable rocks. 
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2. 2. Change in remanent magnetization during storage in the laboratory 

After the first measurement on the natural remanent magnetization was made, all 
the specimen cores were left in the laboratory in such a way that the present geomag- 
netic field keeps a definite angle with the first observed direction of remanent magneti- 


zation. Re-measurement of the remanent magnetization was carried out 30 days after 
the first measurement. 


As for the rock samples having scattered polarization, considerable change in the 
direction of remanent magnetization was found during the storage for times as short 
as a month. The results illustrated in Figs. 2—a, -b, -c are those for the samples 2502, 
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2503 and 2601 respectively. As is seen in these figures, the changes took place in such 
4 manner that direction of remanent magnetization tends to approach that of the applied 
geomagnetic field along the great circle paths. It must be also noticed here that the 
amounts of the change vary core by core but in general they are smaller for the 
specimen cores with reverse polarization (up-dip) which are shown by the hollow circles 
in Fig. 2. As to the sample 2401 with well-concentrated polarizations we could not 
detect any appreciable change in both the direction and the intensity of the remanent 
magnetization during storage in the laboratory for six months. 
2. 3. Demagnetization with alternating magnetic field 

An alternating magnetic field was applied to the specimen cores with the aid of a 
solenoid coil which was tuned to the 50 c.p.s. mains by means of a condenser system. 
This coil was placed horizontally in a field free space produced by a pair of Helmholtz 
coils which compensate the horizontal and vertical components of the geomagnetic 
field. The specimen cores were further rotated around vertical axis during the 
demagnetization procedure. Change in the direction and the intensity of the remanent 
magnetization with the successive increase of the peak intensity of the alternating magne- 
tic field was measured. The results are shown in Fig. 3 and Fig. 4. Figs. 3-a, —b, -c 
and -d represent the variation of intensity of magnetization with applied A.C. field for 


the samples 2502, 2503, 2601 and 2401 respectively. In these figures, the relative intensity 


M/M, was adopted for the ordinate, M, being the initial intensity. It is clearly found 
that the remarkably great decrease in the intensity took place at the lower value of the 
A.C. field such as 50~100 Oe in the case of the specimen cores which showed a large 
sift of the direction of remanent magnetization during a storage in the laboratory. It 


may also be significant to note here that the intensity of remanent magnetization — 


increases at the earlier stage of the demagnetization for the specimen cores with the 


reverse polarity. Typical examples of such a change are seen in the specimen cores” 


® 
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Fig. 3. Change in the intensity of the remanent magnetization with progressive 
demagnetization by alternating magnetic field. 
@: thermoremanent magnetization 
©: natural remanent magnetization 


a. sample: 2502 
b. “ = 2003 
c . : 2601 
d i : 2401 


2502: 3-1, 4-2, 2503: 1-1, 4-1. 

Demagnetization of the thermoremanent magnetization produced in the laboratory 
for the same specimens is also shown in Fig. 3 and compared with that of the natural 
remanent magnetization. A distinct discrepancy is seen in the mode of demagnetization 
between these two remanent magnetizations. Since some preliminary experiment showed 
that a heating to about 600°C does not affect the magnetic properties of the ferro- 
magnetic minerals contained in these rocks, this discrepancy is considered to be genuine 
and suggest that the unstable remanent magnetization of these samples zm situ should 
not be of the thermoremanent origin. 

As to the stable rock sample (2401) the demagnetization curves for both N.R.M. 
(natural remanent magnetization) and T.R.M. (thermoremanent magnetization) which 
were illustrated in Fig. 3-d are extremely similar to each other, about 90 % of the 
intensity of remanent magnetization still remaining after the demagnetization by the 
A.C. field of 400 Oe peak value. There may be no room for doubting that the natural 
remanent magnetization of this rock sample is of the thermoremanent origin. 

Considerable amounts of change in the direction of the remanent magnetization 
also took place for the unstable specimen cores through the A.C. field demagnetization. 
Typical examples are shown in Figs. 4-a, -b and -c which correspond respectively to 
the rock samples 2502, 2503 and 2601. It is remarkable that the change occurs usually 
along the great circle path and for some specimen cores the polarization reversed 
through the progressive demagnetization (see the specimen core 2502: 5-1 in Fig. 4-a). 
In the case of the stable rock sample (2401), it may be safely said from Fig. 4-d that 

| the direction of the remanent magnetization is not affected by the A.C. field of which 
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the highest is 400 Oe. 
2. 4. Demagnetization by heating 

Change in both the direction and the intensity of the remanent magnetization by 
heating experiment in a field free space was examined for some specimen cores. 
Typical thermal demagnetization curve for the rock sample 2503 is shown in Fig. 5. 
It may be noticed in the figure that, as 
to the specimen core with reverse polari- M 
ty (2503, 1-2), the relative intensity 
M/M, increases temporarily in the earlier 
stage of demagnetization and then de- 
creases, while as to the specimen core 
with normal polarity (2503, 6-2), it 
decreases somewhat abruptly in the ae 
earlier stage of the demagnetization, 
the intensity being about one half of 
the initial one after the thermal de 
magnetization of 100°C. 

It was also confirmed that the ; 6 SapnliGao EETSsUME ann ee 


direction of the remanent magnetization Temperature 


of these unstable specimen cores changes Fig. 5. Change in the intensity of the remanent 
Aaall 1 - sel h magnetization with thermal demagnetization. 
gradually along the great circle pat sample : 2503 


through the thermal demagnetization 
procedure, especially the specimen core of 2503, 6-2 reversing its polarity at the de- 
magnetization of 190°C. 

These features of the unstable rock samples for the thermal demagnetization 
resemble very much to those for the A.C. field demagnetization. 
2.5. Intensity ratio of NRM against TRM, J dts 

The intensity values of the natural and thermo remanent magnetizations for some 
specimen cores are shown in Table I together with the ratio between them. The 
intensity of NRM of the unstable dolerite sheet is generally of the order of magnitude 
of 107° emu/gr, which is only some several percent of the TRM value. This makes a 
distinction with the stable rock sample (2401), in which the value of J,//T¢ is 0.34. This 
extremely small value of the Jn/JT- ratio in the unstable rock specimens may be the 
most positive evidence showing that the NRM of these specimens is not directly 
attributable to the TRM of the same specimens, provided that the intensity of the 


geomagnetic field has remained unchanged. 


3. Microscopic Observation of the Samples 


c observation was made on both thin sections and polished sections of 
Direct visual examination of the polished sections by reflected 
the general features of 


Microscopi 


the selected rock samples. 
light under an ore microscope yielded valuable information on 
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Table I. Intensity of natural remanent magnetization and thermoremanent 


magnetization of dolerite sheets and andesite lava. 


Specimen no. Ils ITe In/I]Te 
Kingosi dolerite sheet 
2502 1-2 2.0610 . 88.26x10°5 
e.m.u./gr. e.m.u./gr. 0.023 
2502 3-2 4.64 187529 0.034 
2502 4-1 7.20 92.27 0.078 
2502 5-2 4.49 107.91 0.042 
2502 6-2 10.33 104.07 0.099 
2503 2-2 Seog 96.11 0.037 
2503 4-2 6.94 86.70 0.080 


Warabidaira dolerite 


sheet 

2601 3-1 4.41 ELEN) 0.057 
Tengudana_ andesite 

lava 

2401 7-1 80.37 239-5 0.336 


the opaque minerals contained in the magnetically unstable rocks. So far as the 
present specimens are concerned, the texture of the ferromagnetic minerals in the 
unstable rocks is distinguished by the following two kinds of features from those in 
the stable rocks. The one is the presence of a part with a bluish grey colour in the 
reflected light under oil immersion. This colour has been often reported as characteristic 
to maghemite (Ramdohr, 1950). The typical examples of this feature are shown in 
Plate I and II where an interchange of an original magnetite (more correctly speaking, 
titanomagnetite) and maghemite (more correctly speaking, titanomaghemite) is clearly 


seen. The other is an irregular crack over a whole polished sections of the ferro- — 
magnetic minerals, being demonstrated in Plate III. It must be also noted here that: 


oe es features often coexist in the same ferromagnetic minerals. In Plate IV a 
rogragl ob he peli of, a hee Re eee 


Table II. Petrographic character of dolerite sheets of Kingosi and Warabidaira and of andesite lava of Tengudana. 


Specim : Chief constituent Quantity of | Quantity of |Average size|A ize, Exsolution i 
p en no Rock mnterals magnetite** | ilmenite of nEREe of Hinenite lamellae of teeta in Beavite Shaw 
(vol. %) (vol. %) (mm) Gam) a in ae eneeted light) (ct mh.* (crack) 
vol. % 
2502-3. Quartz-bearing | augite, olivine, magnetite, 
Kingosi, Sidara, | olivine-augite- | ilmenite, plagioclase, 3.4 0.5 0.15 0.1 
Aiti pref. dolerite interstitial alkali feldspar, ; ; ; 2.5 gS UES bd 
quartz, (chlorite) 
2502-4 Seba amet augite, olivine (altered) 
olivine-augite- magnetite, ilmenite, 4.2 On Bx 
dolerite plagioclase, alkali feldspar, ; vale 0.1 2.5 rare~common, | rare 
(chlorite), quartz 
2502-6 Quartz-bearing augite, olivine (altered) 
olivine-augite- | magnetite, ilmenite, 3.3 0.5 0.15 0.1 4.5 common rare~common 
dolerite alkali feldspar, quartz, ; ; (mh) 
(chlorite) 
2503-1 Fine-grained augite, magnetite, ilmenite, 
Kingosi, Sidara, quartz-augite- plagioclase, interstitial alkali 1.3 0.9 0.05 0.1 18 absent absent 
Aiti pref. dolerite feldspar, quartz, (chlorite) 
2503-2 Quartz-bearing | olivine (altered), augite, 
olivine-augite- | magnetite, ilmenite 3.4 0.8 0.2 Oil fi common rare 
dolerite plagioclase, quartz, (chlorite) 
2503-3 Fine-grained augite, magnetite, ilmenite, 
quartz-augite- plagioclase Ns? Nar 0.1 0.1 6 rare~common | common 
dolerite (mt, mh) 
2503-4 Quartz-bearing | augite, olivine (altered), 
olivine-augite- | magnetite, ilmenite, Pek 0.5 OZ 0.15 WW absent common (mt) 
dolerite plagioclase, interstitial alkali 
feldspar, quartz, (chlorite) 
2503-5 Trachytic augite, magnetite, ilmenite, 
dolerite (with | plagioclase, quartz, 1.9 Ong 0.1 0.1 8 rare~common | rare 
granophyric interstitial alkali feldspar, 
vein) (chlorite) 
2503-6 Quartz-bearing | augite, magnetite, ilmenite, 
augite-dolerite plagioclase, interstitial 1.6 1.0 OFS 0.15 6 rare~common | rare 
alkali feldspar, quartz, 
(chlorite), (carbonate) 
2601-3 Olivine-augite- | olivine (altered), augite, ie 0.1 5 ee a 
Warabidaira, dolerite magnetite, ilmenite, : 1.8 0. : ; 
the valley of plagioclase, interstitial alkali 
the river Huzi, feldspar, zeolite, (chlorite) 
Yamanasi pref. 
2401-7 Olivine-bearing| phenoc.: augite, hypersthene, 0.01 pireent absent Aheene 
Tengudana hypersthene- She ugh n.d. n.d. 0.01 
i ‘ti pref. | augite-andesite plagioclase _ 
eidara, Aitt Dr g groundm.:augite, plagioclase, 
magnetite, ilmenite, 
es d 2503-4 
- : ; ; ite except 2503-1 an — 
* i]: ilmenite mt: magnetite mh: maghemite ** including maghem1 Pp 
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Plate I. Specimen: 2503-6 (x1120) Plate II. Specimen: 2503-6 (x 1050) 


Titanomagnetite (grey) containing Titanomagnetite (grey) containing 
titanomaghemite (right side, light titanomaghemite (centre, light grey 
grey part). part). 


Plate III. Specimen: 2503-4 (1190) Plate IV. Specimen: 2401-7 (x1820) 
Titanomagnetite (light grey) Fresh titanomagnetite. 
containing ilmenite lamellae 
(right side, grey part) and cracks. 


4. Identification of Titanomaghemite 


Selected rock specimens having unstable magnetization were crashed by a non- 
magnetic gun-metal crasher and the ferromagnetic fraction was separated magnetically 
in a usual way (Akimoto and Katsura, 1959) from silicate minerals. Identification of 
phases present in the ferromagnetic fraction was carried out by an X-ray diffractometric 
method. X-ray diffraction patterns of the separated ferromagnetic fraction were recorded 
by an X-ray diffractometer “Norelco” with Fe K radiation in a fairly large scanning 
speed for 20 angle range from 90° to 20°. From this general survey we could estimate 


roughly the quantity of ilmenite-hematite series minerals which were occasionally 
involved in the specimens in the form of exsolution lamellae. The lattice parameter 
of the spinel phase was determined by the precise measurement of the diffraction angles 
of the specified crystal planes. Practically we measured the diffraction angles 26 of 
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(220), (311), (400), (511; 333), (440) planes by the “ Norelco ” with Fe Ka radiation with 


the scanning speed of : ° per minute. We adopted a mean value of the ones calculated 


(311) (311) 


<—__ titanomaghemite 
<— titanomaghemite 


(tence poner, (ee ae 
20=44.5° 45.5° 20=44.5° 45.5° 
Specimen no. specimen no. 

2503-2 2601-3 


Fig. 6. Example of the Norelco diffrac- 
tion patterns indicating the presence of 
the titanomaghemite. (Fe Ka radiation) 


from these five 20 angles as the lattice 
parameter of the specimen. 

Detailed examination of the diffraction 
chart of the specimens with the extremely 
unstable magnetization (2502-4, —-5, -6, 2503-2 
and 2601-3) revealed distinctly the existence 
of the spinel phase with smaller unit cell 
dimension in addition to the principal spinel 
phase. As was illustrated in Fig. 6, the 
quantity of this secondary phase is less 
than the principal phase. Considering that 
this situation is in a good agreement with 
the experimental evidence that the phase 
with bluish grey colour in the reflected light 
was observed in the microscopic examin- 
ation of the polished section of the con- 
cerned specimens, this secondary phase may 
be identified as titanomaghemite. 

The results of these X-ray analyses 
were summarized in Table III, where the 
values of unit cell dimension and _ the 
estimated quantity of the titanomaghemite 
relative to the principal spinel phase were 
given together with some remarks on the 
general profile of the X-ray diffractometer 
measurement. 

In order to obtain the direct evidence 
of the presence of titanomaghemite, chemical 
analysis was carried out by T. Katsura of 
Tokyo Institute of Technology for three 
main components FeO, Fe,O; and TiO, of 
the ferromagnetic fractions. Five  speci- 
mens: 2502-2, 2502-4, 2503-2, 2503-4 and 
2601-3 were analysed. The results are given 
in Table IV and represented on the diagram 
of FeO-Fe,O;-TiO, ternary system shown 
in Fig. 7. The composition of all the pre- 
sent specimens deviates from TiFe,O,-Fe,0, 
solid solution line and especially that of the 
specimens with extremely unstable magneti- 
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Table II]. Results of the X-ray analyses of the ferromagnetic minerals separated from the 
magnetically unstable dolerite sheet. 


E Me f = 
| pistes cae eee ba. »| Estimated quantity 
Specimen no. a, (principal ay (titano- of titanomaghemite |II-Ht series minerals 
| spinel phase) maghemite) a, /ay 
| Kingosi | ( ee 
dolerite sheet | 
2502-1 | 98.393 A n.d. 0.3 hardly detectable 
2502-2 | 8.304 (Ses 0.6 . 
2502-3 8.401 8.371 0.3 detectable 
2502-4 8.406 8.375 0.2 39 
2502-5 8.408 8.375 0.3 hardly detectable 
| 2502-6 8.405 8.371 0.3 3 
| 2503-2 8.419 8.386 OR2 a 
| 2503-4 8.407 n.d. 0.1 detectable 
| Warabidaira 
dolerite sheet | 
2601-3 | 8.419 8.38 0.2 hardly detectable 


Table IV. Chemical composition, Curie temperature and coercive force of the ferromagnetic 
minerals separated from the magnetically unstable dolerite sheet. 


Chemical composition : ; 
caer Se Se Gate Curie tem-| Coercive 
FeO Fe,0; TiO, +TiO, FeO Fe,0; TiO,/perature force 
inwt. % is ok. % ? in mol. % 
Kingosi 
dolerite sheet 
2502-4 20.01 44.89 7.70 72.60 42.45 42.85 14.70| 565°C 41 Oe 
2502-6 19.41 45.63 18.53 83.57 34.29 36.27 29.44| 580 39 
2503-2 19.89 42.76 20.95 83.60 34.31 33.19 32.50] ( 225 43 
2503-4 16.02 42.42 16.98 15 31.81 37.88 30.31] 530 35 
Warabidaira 
dolerite sheet 
2601-3 16.76 41.22 16.78 74.76 33.26 36.80 29.94} 550 31 


zation is situated near TiFeO;-Fe,O; line. Since the amount of the ilmenite-hematite 
series minerals in the form of exsolution lamellae was estimated from the microscopic 
observation (see Table II) to be less than about ten percent, this composition may be 
approximately regarded as that of the spinel phase with marked vacancy in the metal 
ion site of the crystal structure. It must be also noted here that the spinel phase with 
bluish grey colour in the reflected light should possess the more oxidized composition 
which lie probably in TiO,-TiFeO;-Fe,O, compositional field. This may be easily 
understood if we recollect the observed fact that the texture of spinel phase is usually 
very. complicated under the ore microscope, being composed of not a homogeneous 


phase but generally of more than two phases. 
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TiO2 


Feo : Fe30, Fe203 
Fig. 7. Chemical composition of the ferromagnetic minerals 
contained in the unstable dolerite sheets. 4 


From these considerations on the results of X-ray and chemical analysis, it may 
be concluded that the titanomaghemite is exclusively found in the rock specimen with 
extremely unstable magnetization. The term of titanomaghemite is used in the present 
paper for the Ti-bearing spinel phase having chemical composition near TiFeO,—Fe,O; 
line or in the TiO,-TiFeO;-Fe,0; compositional field in the FeO-Fe,0;-TiO, diagram. 
Recently Basta (1959) used the term of titanomaghemite for some homogeneous spinel 
phase in South Africa having chemical compositions in a triangular field TiO.-TiFeO, 


-Fe,03. ; ~ 7 
5. Magnetic Properties of Titanomaghemite : Saas ? . - 
5. 1. Thermomagnetic analysis Sa, ec te bade 


investigatio ' 


Natural Occurrence of Titanomaghemite and Its Relevance to the Unstable Magnetization of Rocks 107 


maghemite (spinel structure) to il- 
menite-hematite series minerals (rhom- 
bohedral structure) does not take 
place at the temperature range below 
600°C. 


The Curie temperature determined 


0.5 F 


from the thermomagnetic curve is 
given in Table IV. As seen in the Digits 
table the values of the Curie tem- (in Hex = 900 Qe) 


perature of the present specimens are 


confined in the fairly narrow tem- 


——— —E 


0 100 200 300 400 500 600 °C 
Temperature 


perature range from about 520°C to 
580°C. These values do not seem to 
be incompatible with the equal Curie 1.0 
temperature diagram of the _ spinel 
phase in the FeO-Fe,O;-TiO, ternary 
system which was presented in a 
previous paper (Akimoto, Katsura and 
Yoshida, 1957), although the present ube 


specimens possess the very abnormal 


2503-2 
chemical composition for the spinel (in Hex=1200 Oe ) 
phase. 
5. 2. Coercive force 
The precise knowledge of the 0 100 200 300 400° 500 600 °C 
j Temperature 
coercive force of the ferromagnetic 


Fig. 8. Thermomagnetic curve of the 
fractions seems to be very important titanomaghemite. 

in examining the origin of the unstable 

remanent magnetization. Measurement of coercive force was practised by means of 
the magnetic balance for some specimens separated from the rocks with unstable 
magnetization. The detailed method of determining the accurate value of the coercive 
force was described before (Akimoto, 1957). The results of the measurement are also 
listed in Table IV, the values of the coercive force of the present titanomaghemite 
specimens varying from 31 Oe to 43 Oe. It may be of some significance that the values 
found are considerably larger compared with the intensity of the geomagnetic field. 
These values do not seem to show any abnormal character for the natural ferromagnetic 
mineral grains. They are the very reasonable ones considering that the average size 
of each grain of the ferromagnetic fractions is in the order of 100~200 (Akimoto, 1955). 
From these considerations it may be safely concluded that the coercive force of the 


titanomaghemite does not have an essential influence on the unstable magnetization 
of the rocks containing them. | 
6. Interpretation of the Experimental Results 


baths section the probable origin of the unstable magnetization of the dolerite 
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sheets will be inferred from the experimental results given in the preceding sections. 

As was stated in section 2.3 and 2.5, if we recollect the remarkable difference in 
the A.C. field demagnetization between N.R. M. and T.R.M. of the unstable dolerite 
sheets and the small value of J,//r, ratio, it may be obvious that N.R.M. of these 
rocks is not attributable to T.R.M. of the ferromagnetic minerals actually contained 
in them. The mode of the A.C. field demagnetization curve of the unstable specimen 
cores such as 2502: 6-1, 2503: 2-1, 6-1 and 2601: 3-2 may be well interpreted as a 
characteristic mode for the isothermal remanent magnetization (I.R.M.) (see Figs. 3-a, 
-b and -c). The extremely small value of N.R.M. of these specimen cores may also 
be explained as a character of I. R. M. 

The microscopic observation as well as the X-ray and chemical analysis of the 
ferromagnetic minerals brought out several important facts concerning the magnetic 
behaviour of these unstable dolerite sheets. Firstly we could identify the titanomaghemite 
phase exclusively in the unstable rock specimen. The titanomaghemite being produced 
by a severe oxidation of TiFe,O,-Fe;0, series minerals is generally regarded as a 
titanomagnetite having a large number of vacant cation sites in the crystal structure. 
Therefore, the presence of the titanomaghemite seems to be a strong evidence showing 
that the mother rocks have been subjected to a severe alteration by hydrothermal 
solutions or weathering. 

Through such the alteration or weathering, the remanent magnetization of the 
.mother rock should also be modified. It seems likely that the T.R.M. of the titano- 
magnetite probably having been contained in the original rock is destroyed according 
as the chemical reaction from titanomagnetite to titanomaghemite proceeds. When 
the original titanomagnetite was completely altered, the remanent magnetization of the 
rocks observable in the laboratory is only due to the isothermal origin produced by 
the magnetic after effect in the geomagnetic field. It is likely that the instability of 
the remanent magnetization of the dolerite sheets originates from this isothermal 
remanent magnetization. Since the alteration or weathering take place generally at 
temperatures far below the Curie temperature of the titanomaghemite, the production 
of the T.R.M. of the titanomaghemite, can be ignored. 

If such is the case, the details of the experimental results for the unstable magneti- 
zation of the dolerite sheets can be interpreted consistently as follows. The scattered 
polarization of N. R. M. is chiefly attributable to the circumstances that the rock specimen 
was left in a laboratory at a random orientation to the present geomagnetic field before 
the first measurement of the remanent magnetization was made. Considering that 
within only a few weeks the remanent magnetization changes its direction towards 
the ambient geomagnetic field during the storage, the time duration amounting to a 
month or so which elapsed before the first measurement since the time of sampling 
seems to be sufficiently long to have scattered the direction of the remanent magneti- 
zation. 

Some complicated behaviour of the A.C. field demagnetization for the natural 
remanent magnetization is also elucidated by considering a difference in the degree of 


_ 
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the alteration or weathering. As to the rock specimens of which original titanomagnetite 
partly altered to titanomaghemite, the remanent magnetization measurable at present 
in the laboratory is expected to be a superposition of the original T.R.M. and the 
I.R.M. of the titanomaghemite phase. If the original T.R.M. possesses the reverse 
polarity to that of the I.R.M. the intensity of remanent magnetization will increase at the 
earlier stage of the A.C. field demagnetization, and occasionally the polarization of the 
remanent magnetization will reverse through the progressive demagnetization. This 
may be the case for such specimen cores as 2502: 3-1, 4-2, 5-1, 2503: 1-1, 4-1 (see 
sections 2:32): 


fod 


7. Conclusion and Further Problem 


From the experimental results described in this paper it is clear that the scattered 
polarization and the unstable magnetization are closely related with each other. The 
physical origin of these phenomena is chiefly attributable to the circumstances that 
the T.R.M. has been destroyed more or less in the rocks by the chemical reactions 
such as alteration or weathering and consequently the remanent magnetization measur- 
able at present is mostly due to the isothermal origin. 

The greater part of the magnetically unstable igneous rocks so far reported may 
be interpreted in terms of the present mechanism. Hence it is expected that the 
titanomaghemite phase may be found in the unstable rocks such as Karroo dolerite 
(Graham and Hales, 1957, Nairn, 1958) or South American dolerite (Creer, 1958). This 
mechanism may also be applicable to the case of the unstable sediments such as 
reported by Kawai (1954). 

If the unstable part of the remanent magnetization is completely removed by the 
demagnetization with alternating field, heating or a combination of both, it is possible 
to use the unstable rock samples for a palaeomagnetic purpose. But it was ascertained 
in the present study that it is practically difficult to distinguish the stable fossil part 
of the remanent magnetization for the rock samples suffered from a severe alteration 
or weathering. Therefore, it is advisable not to use for the palaeomagnetic purpose 
the rock samples in which titanomaghemite phase can easily be identified by an ore 
microscope. 

The experimental results of this study may have some bearing on the production 
mechanism of chemical remanent magnetization (C.R: M.). It has recently been verified 
experimentally that the chemical processes can generate a stable remanent magnetization 
under the influence of magnetic field even at a temperature much below the Curie 
point (Haigh, 1958, Nagata and Kobayashi, 1958, Kobayashi, 1959). Such stable chemical 
gnetization, however, could not be found in the present altered dolerite 


The fact that in certain cases the chemical remanent magnetization does not 


remanent ma 
sheet. 
grow to a stable and intense magnetization was recently discussed by Kobayashi 
(1960) theoretically. - ee: 
Another important result of the present study is the discovery of the titano- 
maghemite of which chemical composition is located nearly on the TiFeO;-Fe,O; line 
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in the FeO-Fe,O,-TiO, ternary system. Although the occurrence of titanomaghemite 
in igneous rocks has been paid little attention so far, it may be found commonly in 
igneous rocks especially in the sheet-like intrusive body. The hydrothermal solution 
is generally active in the sheet-like instrusive body during and immediately after the 
solidification of the magma. Although some of its magnetic properties: i.e. Curie 
temperature, coercive force, etc., were obtained in the present study, there is still much 
room for understanding its various peculiar magnetic behaviours. One of the problems 
to be solved may be the possible difference in magnetic viscosity between titanomagne- 
tite and titanomaghemite. If it is ascertained that titanomaghemite shows a larger 
magnetic viscosity, the remarkable situation that the remanent magnetization of the 
rocks containing the titanomaghemite changes its direction very easily under the am- 
bient magnetic field would be understood. Systematic study along this line is now 
undertaken. 
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